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Abstract 

We present an experimental study of the isotropic-nematic phase transition in an aqueous mixture 
of charged semi-flexible rods (/c? virus) and neutral polymer (Dextran). A complete phase diagram 
is measured as a function of ionic strength and polymer molecular weight. At high ionic strength 
we find that adding polymer widens the isotropic-nematic coexistence region with polymers prefer- 
entially partitioning into the isotropic phase, while at low ionic strength the added polymer has no 
effect on the phase transition. The nematic order parameter is determined from birefringence mea- 
surements and is found to be independent of polymer concentration (or equivalently the strength 
of attraction). The experimental results are compared with the existing theoretical predictions for 
the isotropic-nematic transition in rods with attractive interactions. 
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I. INTRODUCTION 



One of the fundamental notions of the theory of hquids, dating back to van der Waals, is 
that hquid structure is determined by the repulsive part of the intermolecular potential. The 
attractive part of the potential determines the density of a liquid by providing a cohesive 
background energy that is largely independent of a particular configuration of molecules 0,0|. 
This is true as long as the liquid is far from its critical point. Due to this reason there has 
been a substantial effort over the past 50 years to use hard spheres as a reference system 
to understand the behavior of all simple liquids^]. Parallel to these endeavors, a theory 
of a liquid of rods with purely repulsive anisotropic interactions was developed by Onsager 
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It was shown that this system exhibits an isotropic-nematic (I-N) phase transition. 
The Onsager theory is based on the realization that the virial expansion of the free energy 
converges for hard rods with sufficiently large aspect ratio at the I-N phase transition, 
in contrast to spheres where the virial expansion fails to describe hard spheres at high 
concentration. 

Once the behavior of a hard particle fluid is understood it is possible to study the influence 
of attractions via a thermodynamic perturbation theory 0, y| . For hard spheres this is 
relatively easy due to the fact that attractions provide an structureless cohesive energy. In 
contrast, extension of the highly successful Onsager theory valid for rods with short range 
repulsions to a system of rods with longer range attractive interactions is much more difficult. 
The difficulties stem from the fact that attractive rods are in their lowest energy state when 
they are parallel to each other. These are exactly the configurations that need to be avoided 
if the second virial term on which the Onsager theory is based, is to accurately describe 
the system Q. In one study the Onsager functional has been straightforwardly extended 
to include an additional attractive interaction (from now on called Second Virial Theory 
with Attraction (SVTA)). Because of the problems already mentioned the author argues 
that SVTA is valid only for very weak attractions. Indeed, at high strengths of attraction 
un-physical states such as a collapse to infinitely dense state are predicted. The physical 
picture that emerges from the SVTA theory is that of a van der Waals like liquid of rods 
where its primary structure (i.e. nematic order parameter) is determined purely by the 
repulsive interactions, while attractions serve as a uniform structureless glue holding the 
rods together at a given density. 



For the theory to work at all densities the free energy of the unperturbed liquid of rods 
needs to take into account third and higher virial coefficients. An alternative theory that 
accomplished this uses scaled particle free energy of hard rods as a basis to study the influence 
of attractive interactions on the I-N phase transition y,y](from now on called Scaled Particle 
Theory with Attractions (SPTA)). The scaled particle expression for hard rods includes third 
and higher virial coefficients. Therefore it is reasonable to expect that this theory would be 
more accurate at higher rod and/or polymer concentrations. An additional advantage of the 
SPTA theory is that it does not assume that the depletion interaction is pairwise additive. 
Computer simulations have shown that pairwise additivity of the intermolecular potential 
assumption is not an adequate approximation when the radius of the polymer is larger then 
the radius of the colloid P, 13 • 

In this paper we experimentally study the influence of attractive interactions on the I-N 
transition and compare them to both SVTA and SPTA theory. As a reference system we use 
an aqueous suspension of semi-flexible rod-like fd viruses. Previous work has shown that the 
behavior of fd virus is consistent with the theoretical predictions for semi-fiexible rods with 
purely repulsive interactions [l2 , ^ 14 1 . Strictly speaking, fd forms a cholesteric phase 
and undergoes an isotropic-cholesteric transition, but because the free energy difference 
between a cholesteric and a nematic phase is small we refer to the cholesteric phase as 
nematic in this paper. Additionally, introducing a finite flexibility into hard rods significantly 
alters both the location and nature of the isotropic-nematic and nematic-smectic phase 
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transition [13, ll^- Here we show that flexibility also changes the isotropic-nematic phase 
transition in rods with attractive interactions. 

We induce attraction experimentally by adding a non-adsorbing polymer to the colloidal 
suspensions, which leads to the depletion interaction where the range and the strength of 
the attractive potential is controlled by the polymer size and concentration respectively 3| • 
Although this work specifically deals with a colloid/polymer mixture its results are of a 
general significance to other anisotropic fluids which have attractive interactions due to 
other reasons (i.e van der Waals attractions). The main difference between polymer induced 
depletion attractions and attractions due to van der Waals forces is that in the depletion case 
there is partitioning of the polymer between coexisting phases |17|. Therefore the strength of 
the interaction between two colloids depends on the phase in which the colloids are located. 

There have been previous experiments on the influence of polymer on the I-N phase 



transition in mixtures of boehmite rods and poly styrene polymers and mixtures of cellulose 
nanocrystals and Dextran polymers Q, Q, Q- 

Other studies related to our work have 
focused on the condensation of rod-like polymers due to the presence of polymer and/or 
multivalent cations 21, The conditions in those studies correspond to the upper left 
corner of the phase diagram in Fig. El We also note that at very high polymer concentrations 
the fd system exhibits a direct isotropic-smectic coexistence and a number of metastable 
complex structures associated with this transition have been described elsewhere 2^, 24 1. 



In this paper we limit ourselves to the I-N transition. In Section II we present the ex- 
perimental details of our measurements. In section III the effective intermolecular potential 
acting between two rod-like particles is discussed. In section IV we present the measured 
phase diagrams as a function of ionic strength and polymer size, and in section V we present 
our conclusions. In the appendix we provide the formulas necessary to calculate the phase 
diagrams in the SPTA theory. 



II. MATERIALS AND METHODS 

Bacteriophage fd was grown and purified as described elsewhere Fd is a rod-like 

semifiexible charged polymer of length 0.88 /im, diameter 6.6 nm, persistence length 2.2/im 
and surface charge density of le~/nm at pH 8.15. All samples where dialyzed against 
20 mM Tris buffer at pll=8.15 and NaCl was added until the desired ionic strength was 
achieved. Dextran and FITC-Dextran with molecular weights (MW) of 500,000 and 150,000 
g/mol (Sigma, St. Louis, MO) were used as the non-absorbing polymer and dissolved in 
the same buffer solution. The samples are prepared in the two phase region of the phase 
diagram as is shown in Fig. ^ Concentrations of coexisting phases were measured using 
absorption spectrophotometry. The optical density of fd is ODag^n'm'^ = 3.84 for a path 
length of 1 cm. To determine the concentration of Dextran polymer we used a mixture of 
95% Dextran and 5% FITC labelled Dextran. The optical density of FITC-Dextran was 
determined by dissolving a known amount of polymer in a buffer solution and measuring 
the OD at 495nm. The relationship between the radius of gyration Rn of Dextran and its 
molecular weig ht (MW) in units of g/mol is Rg[A\ = 0.66(MW)°-^3 ^26^. The reason for 
the small exponent 0.43 is due to the fact that Dextran is a branched polymer. The volume 
fraction of polymer 0poiymer was calculated by 0poiymer = pf^-^g; where p is the polymer 



number density. The order parameter of the nematic phase was measured with a Berek 
compensator, by placing the suspension into a quartz x-ray capillary with a diameter 0.7 
mm (Charles Supper, Natick, MA). Samples were aligned with a 2T magnetic field [2^ and 
the birefringence was measured. The order parameter [S) is obtained using the relationship 
An = SpMfio where is the number of rods per unit volume of fd virus. An is the 
birefringence measured using Berek compensator on an Olympus microscope and S is the 
nematic order which varies between for the isotropic phase and 1 for a perfectly aligned 
phase. The birefringence of perfectly aligned fd, hq = 3.8 x 10^^ ± 0.3 x 10~^ ml/mg, was 
recently obtained from x-ray experiments 



III. INTERMOLECULAR POTENTIAL 



When a colloid is suspended in a polymer solution it creates around itself a shell from 
which the center of mass of a polymer is excluded. When two colloids approach each other 
there is an overlap of the excluded volume shells which leads to an imbalance of the osmotic 
pressure that is exerted on each colloid. This results in an effective attractive potential known 
as the depletion potential. In the Asakura-Oosawa model (AO), polymers are assumed to 
behave as spheres of radius R^^, which can freely interpenetrate each other but interact 
with colloids via hard core repulsive interactions [16]. The relationship between R^'^ and 
radius of gyration of a polymer (Rg) is i?^^ = 2Rg/ ^/^^. This approximation is valid as long 
as the size of the colloidal particle is much larger then the radius of the penetrable sphere 
If the size of a colloid is equal to or smaller then R^'~', the colloid can 
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penetrate into the open polymer structure without overlapping any of the polymer segments. 
In this case the range and the depth of the attractive depletion potential will be significantly 
weaker when compared to the predictions of the AO model. In our experiments the diameter 
of the polymer is up to 5 times the diameter of the rod-like virus and therefore we expect that 
the depletion potential significantly deviates from the Asakura-Oosawa penetrable sphere 
model. 

Since there are no analytical results on the depletion potential between rod-like colloids 
we estimated it using computer simulations. The method used to obtain the potential is 
described in detail in the paper by Tuinier et. al. ; here we briefiy outline the procedure. 
Two spheres, cylinders or walls are set at a fixed distance apart and an attempt is made 



to insert a non-self-avoiding polymer molecule at random positions. When simulating the 
depletion potential between the cylinders they are oriented in perpendicular directions. If 
any segment of the polymer overlaps with either colloid, the insertion attempt fails and the 
polymer is not counted. The profile of the depletion potential is then equal to 

f/dcpiction(a;) = kBT{N{oo) - N{x)) (1) 

where N{x) is the number of polymers successfully inserted in the simulation box when two 
colloidal objects are a distance x apart. A^(oo) is the number of polymers inserted when two 
colloids are apart at a distance which is much larger then the range of the intermolecular 
potential. 

The depletion potentials between walls, spheres and rods obtained from the simulations 
are shown in Fig. |21 From the exact results, it is known that the depletion potential at 
small separations between two parallel walls induced by AO penetrable spheres is equivalent 



to the depletion potential induced by polymer (without excluded volume interactions), if 

If we use this fact, the simulation results for the depletion potential 
between two plates (indicated by open circles in Fig. ^ are in a very good agreement with 
the potential predicted by the AO theory (indicated by the full line in Fig. E)), as long as 
the separation between the plates is smaller then 3Rg/2. At larger separations we observe 
that the potential exerted by the polymer has longer range attraction than the equivalent 
penetrable sphere, as was previously noted 31[. This is because a polymer is only spherical 
on average and will adopt elongated conformations on occasion. The simulation results 
for the depletion potential between two spheres immersed in a polymer suspension with 
Rg / RcoWoiA = 3.36 is significantly weaker then what is predicted by the penetrable AO sphere 
model. The reason for this is that a small sphere has a high probability of penetrating a 
polymer with a large radius of gyration since polymers have very open structures. The rods 
have a profile of an infinite plane in one direction and a profile of a sphere in the other 
direction. It follows that a cylinder with the same diameter as a sphere is less likely to 
interpenetrate with a polymer coil. Therefore the depletion interaction between cylinders is 
stronger than between spheres of equal diameter and weaker then the depletion interactions 
between two walls. Even for the case of cylinders, the potential obtained from the AO model 
significantly overestimates the strength of the potential obtained from the simulation as is 
shown in Fig.|21 In this paper we assume that the strength of the depletion potential between 



two cylinders oriented at an angle 7 scales as 1/ sin 7, but the shape remains independent 
of 7. To verify this we have simulated the potential between two cylinders that are either 
parallel or perpendicular to each other. For these two cases we obtain depletion potentials 
that are almost identical to each other after they are rescaled by a constant. This supports 
our assumption that the shape of the depletion potential is independent of the cylinder 
orientation. 

The total interaction potential between two fd viruses in a /c?/Dextran mixture is a 
combination of hard core repulsion, a steep short range electrostatic repulsion, and the 
longer range depletion attraction described above. As the ionic strength decreases both the 
range and the depth of the potential decreases as is shown in Fig. IHl This is due to the fact 
that Dextran is an uncharged polymer and therefore the depletion attraction is independent 
of the ionic strength. Decreasing the ionic strength results in longer range electrostatic 
repulsion which screens out ionic strength independent depletion attraction. 

The short range electrostatic repulsion and longer range depletion attraction scale as 
1/ sin(7) where 7 is the angle between two rods. Therefore, to a first approximation the 
position of the minimum of the intermolecular potential does not change when the angle 
between two rods changes; only the magnitude of the minimum changes. To account for the 
rapidly decaying electrostatic repulsion we re-scale the hard core diameter to an effective 
hard core diameter D^s as was described previously We note however that the use 

of Dcs is rigourously justified only in the dilute regime where the second virial coefficient 
quantitatively describes the system, ie. at the isotropic-nematic transition of pure rod 
suspensions. Therefore one of the causes of the discrepancy between theory and experiments 
stems from our crude treatment of the electrostatic interactions. As discussed previously the 
attractive depletion potential is weaker then the predictions of the AO model. To account 
for this in the calculation of the phase diagram we simulated the depletion potential for 
experimentally relevant parameters. The simulated potential is mapped onto AO model 
where the effective concentrations of the interpenetrable spheres (pes) and effective polymer 
radius (-R^p) are adjusted to yield the best fit to the simulated potential. We define p^s = O-p 
where p = N/V is the actual number density of AO penetrable hard spheres and R^^ = 
PR^^. Surprisingly we find that the a is much smaller then 1 while /3 is only slightly smaller 
than 1 for parameters used in our experiments. This can be seen in Fig. Elwhere the range of 
the depletion potential between two sphero cylinders for the AO penetrable sphere model is 



almost identical to the range of the simulated potential while the depth is very different. The 
reason for this is that the AO model underestimates the depth of the depletion attraction 
at large distances (Fig. 12)). If the polymer size is increased further we observe that the value 
of jS will start decreasing rapidly. The comparison between the simulated potential and the 
effective potential used in the theoretical calculations of the phase diagrams is shown in the 
inset of Fig. El The phase diagrams are calculated using the effective rod diameter Des, the 
effective polymer radius -R^, and the effective polymer concentration pefr- The calculation 
of the SPTA and SVTA phase diagrams for semi-fiexible rods is described in the Appendix. 
Once the phase diagrams are obtained the polymer concentrations are rescaled back to the 
actual volume fraction of polymer. Specifically we calculate the theoretical phase diagrams 
using pcff and then in order to compare with experiment we plot the theoretical results using 
P = Pcff/a- 



IV. RESULTS 

In Fig. ^ a typical phase diagram for a mixture of hard rigid rods and polymers is 
indicated by thick full lines. As was shown in previous work by Lekkerkerker et. al. 
adding polymer widens the isotropic-nematic coexistence and leads to partitioning of the 
polymer between isotropic and nematic phases. In the same figure dashed lines indicate the 
phase diagram of a mixture of semi-fiexible rods and polymers. The influence of the flexibility 
on the isotropic-nematic pha se transition is well studied for the case of rods with hard core 
repulsive interactions 33|,|3^. Flexibility increases the concentration of the I-N co-existence, 
decreases the width of the I-N coexistence, and reduces the order parameter of the nematic 
phase co-existing with the isotropic phase. In Fig. |3] the theoretical phase diagrams for two 
equivalent systems of rods with attractions are shown with the only difference being the 
flexibility of the rod. For the case of the rigid rods the concentration of the polymer needed 
to induce widening of I-N coexisting phases is much lower than for that of semi-flexible 
rods. This is due to the fact that to compress semi-flexible rods, the polymer has to work 
against both rotational and internal bending contributions to the entropy. We conclude that 
flexibility also suppresses the formation of the nematic phase in attractive rods. Next we 
proceed to compare the theoretical phase diagrams to experiments. 

A representative experimental phase diagram of an /c?/Dextran mixture at high ionic 



strength is shown in Fig. El Two features of the phase diagram are in quahtative agreement 
with the theoretically predicted one. First, introducing attractions widens the isotropic- 
nematic coexistence. Second, at intermediate polymer concentrations polymer preferential 
partitions into the isotropic phase. At very high polymer concentration the rods and poly- 
mers are essentially immiscible with a nematic phase of pure rods coexisting with an isotropic 
phase of pure polymers. This part of the phase diagram has been measured elsewhere j2J]. 

We proceed to study the influence of the ionic strength on the phase behavior. Changing 
ionic strength significantly modifies the interaction potential as was shown in Fig. IHl The 
phase diagrams at three different ionic strength are shown in Fig. IHl The experimentally 
measured phase diagram at 50 mM ionic shows that the addition of the polymer has no 
effect on the coexistence concentrations of the I-N transition. This is in stark disagreement 
with theory which predicts strong partitioning of the polymer. The implication from these 
experimental results is that the depletion attraction is completely screened by the long 
range electrostatic repulsion. In calculating the potential energy between charged rods in 
the presence of neutral polymer we are summing two large terms of opposite signs (Fig. Oj). 
Small inaccuracies in the theory of either of these terms could account for the discrepancy 
between the theoretical and experimentally observed phase diagrams. 

As the ionic strength is increased to 100 mM the addition of the polymer initially increases 
the width of the co-existence concentration, while at very high polymer concentration we 
observe re-stabilization of the I-N transition. This was also observed in mixtures of fd and 
Dextran (MW 500,000) at 100 mM. This observation can be explained by the fact that 
restabilization of the I-N phase transitions occurs when the polymer is in the semi-dilute 
regime. In this regime the range of the depletion interaction is of the order of the correlation 
length (polymer mesh size), which is smaller then the radius of gyration [35]. Moreover 
the correlation length decreases with increasing concentration. Since the range of attraction 
decreases in the semi-dilute regime, the long range electrostatic repulsion will screen out any 
depletion attraction in the semi-dilute regime. This mechanism of depletion re-stabilization 
was previously observed in mixtures of charged spherical colloids and polymer mixtures in 
aqueous suspension [36|. 

At the highest ionic strength of 200 mM a relatively low concentration of polymer is 
needed to induce a complete phase separation between a polymer-rich, rod-poor isotropic 
phase and a rod-rich, polymer-poor nematic phase. At this ionic strength no reentrant I-N 



phase behavior is observed for all accessible polymer concentrations. The phase behavior 
at this ionic strength has the same qualitative behavior as predicted by both SVTA and 
SPTA. However, when the theory is quantitatively compared to experiment there are large 
differences between the predicted phase boundaries (full lines) and experimentally measured 
phase boundaries. Perhaps the fact that the disagreement between theory and experiment 
is worst at low ionic strength and high rod concentration indicates that our approximation 
of treating electrostatically repulsive rods as hard rods with an effective diameter D^s is 
invalid under these conditions, as previously mentioned. 

If attractions are introduced to a hard sphere system, the assembly will decrease its 
energy by decreasing the average separation between spheres, which in turn increases the 
density of the stable liquid phase. Unlike spheres, rod-like particles with attraction have a 
more complex interaction potential. They can lower their interaction energy not only by 
decreasing their separation (increasing their density), but also by decreasing their relative 
angle (increasing their order parameter), or a combination of both. To distinguish between 
these possibilities we measure the order parameter of the nematic phase in coexistence with 
the isotropic phase as is shown in Fig. [7| In Fig. [7^ we compare the order parameter of rods 
with and without attractive interactions (ie. with and without the addition of polymer) and 
find that the order parameter is determined by the concentration of rods only. The nematic 
order parameter is plotted as a function of polymer concentration in Fig. [3d to illustrate the 
independence of the coexisting nematic order parameter with increasing attraction. This is 
further confirmed with the graph in Fig. |H1 which shows that the nematic order parameter is 
independent of polymer concentration even well into the nematic phase. This is in agreement 
with both the SVTA and SPTA theories which predict that the measured order parameter 
of the nematic rod/polymer mixture will depend only on the concentration of rods and be 
independent of the level of attraction (ie. polymer concentration). We note that the order 
parameter data is noisy because of the intrinsic high viscosity of the fd/Dextran solutions. 
This high viscosity makes it difficult to create nematic monodomains even in magnetic fields 
up to 8T. 



V. CONCLUSIONS 



We have presented quantitative measurements of the isotropic-nematic phase transition in 
a binary suspension of rod-hke particles (fd) and spherical polymers (Dextran). The widen- 
ing of the coexistence concentrations and partitioning of the polymer predicted theoretically 
are observed in these experiments on /(/-Dextran mixtures. As discussed in the previous 
paragraph, our measurements indicate that the liquid of rods can be thought of as a van 
der Waals liquid where the order parameter of the nematic phase is determined by repulsive 
interactions, while attractive interactions provide structureless cohesive energy. Within the 
admittedly noisy experimental data, we find that the order parameter is determined solely 
by the rod concentration and not by the polymer concentration, or equivalently, the strength 
of attraction. However, even after taking the following effects into account; the possibility 
of the virus and polymer interpenetrating, the charge of the virus, and the semi-flexibili ty o f 
the virus we found large quantitative differences between the theory and the experiment j22| ■ 
Notably, the theory severely overestimated the strength of the polymer induced attraction. 
The difference is especi ally pronounced in the nematic phase and at low ionic strength. This 



and previous work [23|, |38| suggest that much remains to be done before we are able to un- 
derstand and predict the behavior of rods whose interactions are more complex then simple 
Onsager-like hard rods. 
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In this appendix we present t 
mixture using the SPTA theory 



le calculation of the phase diagram for a rod-polymer 
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Several misprints in the original paper are 



The approximate free energy of the colloid-polymer mixture is given by 



the following expression 



17|: 



Fc+p(0) = - Hp < yfree(0) > (2) 

where Fc{4>) is the free energy of colloid suspension at volume fraction 0. The co-existence 



concentrations for the I-N transition predicted by the scal ed p article theory are in very close 
agreement with the results from the computer simulations [41]. This indicates that the scaled 
particles theory provides a good approximation for third and higher virial coefficients. The 
system is assumed to be in equilibrium with a polymer reservoir which is separated from the 
colloid-polymer mixture by a membrane permeable to polymers only. The osmotic pressure 
of the polymers in the reservoir is Up = pksT, with p the polymer number density. \4i.ee 
is the free volume available to a polymer in a solution of pure hard particle colloids. It 
is assumed that Vfree in a polymer /colloid mixture is equal to the Vfree in the pure colloid 
suspension. In this sense Eq. El is a thermodynamic perturbation theory. 

The expression for the free energy of a pure hard sphero cylinder colloidal suspension is 
given by the scaled particle theory developed by Cotter j^] '■ 
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+ -n3(5,a) - 



where is the volume fraction of spherocylinders 
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The coefficients 112 and are given by the following expressions 
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(35 - 1)2 (35 - 1)2 
and the parameter 5 is the overall length to diameter ratio of the spherocylinder 5 = 
The function o"(a, L/ P) is an expression that accounts for the rotational entropy of the rods 
and the entropy associated with the loss of configurations due to confinement of the bending 
modes of the semi-fiexible rods in the nematic phase has been derived by extrapolating 



between the hard rod and the fiexible chain limits 
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45|. In this paper the expression 



for a obtained by Hentschke is used for numerical calculations and is given by 



p) = ln(a) - 1 + 7re " + —(a - 1) + — In (^cosh ^ 



L a — 1 



(7) 



The function C,{a) that describes the interactions between rods at the level of second 
virial coefficient is given by: 



^, , 2/2(0;) 

For this calculation we assume the Onsager ansatz for the orientational distribution func- 
tion given by : 

^/ /^xN acoshfa cos(6')) . . 

/(a..cosW)^ 4jnh(«) ■ 

The expression for the free volume in a sphero cylinder suspension is given by: 



u{<P,6,q) = Y^ = (l-0)exp I - I A{6,q) { ]+B{6,q) [ ) +C 
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where 
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The ratio of the polymer diameter to the rod diameter is given by the parameter q. After 
the expression for the scaled particle free energy ( E)) is obtained, we insert the Onsager 
approximation for the orientational distribution functions /(«) and minimize the free energy 
at different rod concentrations with respect to the parameter a to find the order parameter of 
the nematic phase at that concentration. To find the concentrations of rods in the coexisting 
isotropic and nematic phases we solve the conditions for the equahty of the osmotic pressure 
and chemical potential. The expressions for the osmotic pressure and the chemical potential 
are : 



n 



,9z/(0)' 



fl = Fc[<p) + h UpX 



(12) 



where is the polymer volume fraction, ^nR^p, and A is the ratio of sphero cylinder volume 
to polymer volume 

A = Afl + ^l^-l)) (13) 



gH 2 

The phase diagram is calculated by first minimizing the free energy with the respect to the 
parameter a and then solving coexistence equationsf fT^ . The SVTA phase diagrams for 
rigid rods are calculated following the calculation of Warren ^7^ . To extend this calculation 
to semi-flexible rods the orientational entropy term in the Onsager free energy was replaced 
by the confinement entropy of semi-fiexible polymers as shown in Eq. [7| To account for 
electrostatics the rod diameter D is replaced with -Deff^|. To correct the AO theory for 
the case of a polymer radius larger than the colloid radius we replace the polymer density 
p with peff and the polymer radius Rg with R^g as described in the text and in Figures 3,5 
and 6. 
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FIG. 1: (Color online) A sequence of images illustrating the preparation of a sample which was used 
in determining the phase diagram, a) A nematic liquid crystal of fd virus in buffer between crossed 
polarizers showing disordered birefringent domains, b) A highly concentrated solution of Dextran 
labelled with yellow fluorescein is added to the transparent fd nematic liquid crystal, c) After 
the sample in (b) is vigorously shaken it phase separates into the coexisting nematic and isotropic 
phases. The macroscopic phase separation takes from few hours to couple of days depending on the 
location in the phase diagram, d) Same sample as image (c) but taken between crossed polarizers. 
Image shows dense birefringent nematic phase on the bottom and Dextran rich isotropic phase on 
the top which is yellow in appearance. Images (a) to (d) are taken on the same sample. 




FIG. 2: Depletion potential U between two walls, two perpendicular eylinders and two spheres 
obtained from computer simulation are shown by open spheres, squares and triangles respectively. 
In the two wall simulation the wall size was 313 x 313 and periodic boundary conditions where 
used. The diameter of the spheres and cylinders is 66 A while Rg of the polymer is 111 A. 
The lines indicate depletion potentials as predicted by the penetrable sphere (AO) model. The 
separation x is the closest distance between two surfaces. The number concentration of the polymer 
p is equal to the overlap concentration p = 3/{4TrRg), while the radius of the penetrable spheres 
is R^^ = IRgl^Jii = 125 A. The AO theory overestimates the potential between spheres and 
between cylinders because polymer deform around colloidal particles. 
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FIG. 3: Total interaction potential C/ as a function of separation x between two viruses {D = 
66 A) oriented at 90° with respect to each other and immersed in a suspension of polymers of 
concentration p = 3/(47ri?^) and radius Rg = 111 A at three different ionic strengths. The 
interaction potential is a sum of electrostatic repulsion and depletion interaction. The effect of 
electrostatic repulsion for fd with net linear charge density 1 e~ /Aof is accounted for by treating 
the the fd as a hard particle with a larger effective diameter L>^g|lll. [3^. Filled circles indicate the 
depletion potential obtained from Monte Carlo simulation of polymers without excluded volume 
interactions. Since the polymer diameter is larger then the rod diameter, the polymer and rod- 
like viruses can easily interpenetrate. This results in a effective depletion attraction which is 
smaller than what is predicted by the AO model (indicated by the full line). The phase diagrams 
corresponding to these interaction potentials are shown in Fig. |SJ Inset : In theoretical calculations 
we approximate the intermolecular potential between rods with an effective hard core diameter 
-CefflUl and attractive potential. The attractive part of the potential is modeled by AO penetrable 
spheres whose effective radius and concentration best fits the potential obtained through computer 
simulation. This effective intermolecular potential is compared to the potential obtained through 
the computer simulation in the inset. In the inset /Oeff/p = 0.39 and R^^ /R^'^ = 0.99. 
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FIG. 4: Phase diagram for rigid and semi- flexible rods calculated using the SPTA theory. The 
boundary between the isotropic(I)-nematic(N) two phase region and the region where a single phase 
is stable is indicated by the thick dashed line for semi-flexible rods and thick full lines for rigid 
rods. Tie lines between the coexisting phases are shown by thin lines. For the flexible particle the 
ratio of the contour length to persistence length is L/P = 0.4. The phase diagram was calculated 
using (5 = 84 and q = 2.2. The polymer concentration is defined as follows copolymer = p—^- 




c^[mg/ml] 

FIG. 5: Phase diagram for a mixture of fd virus and Dextran (MW 500,000, Rg = 176 A or 

ionic strength. The measured points indicate the rod and polymer 
concentrations of the coexisting isotropic and nematic phases. The full line is a guide to the eye 
indicating the two phase region. Tie-lines are indicated by thin full lines. The SPTA and SVTA 
predictions are indicated by the dotted lines and dashed lines respectively. 
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FIG. 6: Phase diagrams of a mixture of fd virus and Dextran polymer (MW 150,000, i?g = 111 A 
or R^'^ = 125 A) at 50 mM, 100 mM and 200 mM ionic strength. Coexisting phases are indicated 

by open circles while the full line is an eye guide separating two phase region from isotropic and 
nematic phases. The predictions of the SPTA and SVTA theories are indicated with dashed and 



dotted lines respectively. The polymer concentration is defined as follows <;!!>poiymer = P 
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FIG. 7: Measurements of the order parameter of the nematic phase in coexistence with the isotropic 
phase for a mixture of fd (rod) and Dextran (M.W. 500,000, polymer) at 100 mM ionic strength. 
Order parameter (S) is graphed as a function of /d concentration (a) and Dextran concentration (b) 
for the coexisting nematic concentrations shown in (c). The order parameter is double valued in (b) 
because along the nematic branch of the coexistence curve there are two different rod concentrations 
with the same polymer concentration as shown in (c). Data in (c) is the same as that shown in 
Fig. [5J Dashed line in (a) indicates the theoretical dependence of the order parameter on the 
concentration of rods as obtained using scaled particle theory. This relationship agrees well with 
experimental data for fd at high ionic strength using x-ray scattering 2^. Arrows in (b) and 
(c) indicate the direction of increasing fd concentration. Below c/'poiymer ~ 0.2 the nematic fd 
concentration is essentially constant. 



1.0 



■ fd concentration ~ 33mg/ml 



0.9 



S 



OS 



0.7- 




0.0 



0.1 




02 



03 



FIG. 8: The order parameter of the nematic phase of the /d/Dextran (M.W. 500 000) mixture 
at 33 mg/ml fd and 100 mM ionic strength as a function of increasing polymer concentration. 
The horizontal line drawn is a guide to the eye showing the independence of the nematic order 
parameter with polymer concentration. The vertical line indicates the location of the nematic- 
isotropic transition. 



